pro realhl2a, var=var,ntime=ntime,output=output

path="data"
g=file_import("data/circular.grd.hl2a22535.831ms.68x256.0.65to1.1.nc")
data=collect(path=path,var=var)
IF NOT KEYWORD_SET(ntime) THEN ntime = 0
IF NOT KEYWORD_SET(output) THEN output='false'
f=data[*,*,0,ntime]
f=reform(f)
xrag0=0.65 ;1.935
xrag1=1.1 ;2.065
yrag0=0.0
yrag1=0.5
zrag0=min(f)
zrag1=max(f)
chars1=3
thick1=2
thick2=2

casename=''
titlename=''
xtitlename="R"
ytitlename="Z"

safe_colors,/first
!x.margin=[10,2]        ;;
!y.margin=[8,20]       ;; 
loadct, 39             ;;
;loadct, 30
!p.background=-1       ;;
!p.color=1             ;;
window,0,xsize=900,ysize=700,xpos=0

if output eq 'true' then begin 
SET_PLOT,'PS'
DEVICE, file=filename,/color,/landscape
endif
contour2,f,g.rxy,g.zxy,/fill,nlevel=10,charthick=thick2,chars=chars1,thick=thick1,$
title=titlename,ytitle=ytitlename,xtitle=xtitlename,$
xthick=thick2,ythick=thick2
if output eq 'true' then begin
DEVICE, /close
SET_PLOT, 'X'
endif


end
